Robustness of the adiabatic quantum search 
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The robustness of the local adiabatic quantum search to decoherence in the instantaneous eigen- 
basis of the search Hamiltonian is examined. We demonstrate that the asymptotic time-complexity 
of the ideal closed case is preserved, as long as the Hamiltonian dynamics is present. In the special 
case of pure decoherence where the environment monitors the search Hamiltonian, it is shown that 
the local adiabatic quantum search performs as the classical search. 
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Although the adiabatic approach to quantum compu- 
tation P, I3] seems to differ significantly from the tradi- 
tional circuit model, it has been proved that these two 
models are, in a certain sense, equivalent Q. However, 
this equivalence does not concern the robustness to noise, 
relaxation, or decoherence. Since the adiabatic schemes 
operate close to the energy ground state it seems natural 
to guess that the adiabatic quantum computer should be 
robust against relaxation effects 2J. The alleged resis- 
tance to noise has been examined by analytic means in 
Ref. 4] and it has been argued that adiabatic quantum 
computers should be robust to decoherence 0, la] ■ Uni- 
tary control errors and resistance to decoherence have 
been numerically investigated in Ref. 6] . 

In this paper, we examine the local adiabatic search 
algorithm Q in the presence of decoherence in the 
instantaneous energy eigenbasis ^ . We demonstrate an- 
alytically a robustness to this particular form of decoher- 
ence in the sense that the asymptotic time-complexity of 
the ideal closed case is preserved, no matter how small 
the Hamiltonian contribution is to the dynamics. Only 
in the wide-open case (lOj . where the Hamiltonian part 
is completely absent, there is a difference in the time- 
complexity. 

Adiabatic quantum computation works by keeping the 
system close to the ground state of a time-dependent 
Hamiltonian. This feature is in contrast with, e.g., holo- 
nomic implementations of quantum gates TTj, which 
share the feature of adiabatic evolution, but where it is 
essential that the gate can operate on arbitrary superpo- 
sitions without too large errors. For the functioning of 
the adiabatic quantum computer in the presence of deco- 
herence, on the other hand, it is sufficient to require that 
the probability of finding the system in the instantaneous 
ground state of H{s) is conserved. This can be seen as 
one possible generalization of the concept of adiabaticity 
to open systems. In this generalized sense the wide-open 
case has an adiabatic limit, although the Hamiltonian dy- 
namics is absent. One may note that the wide-open case 
can be seen as a quantum computational scheme in its 



own right, a "wide-open adiabatic quantum computer", 
where the dynamics is governed by pure decoherence. A 
different approach to the concept of adiabaticity for open 
system has been put forward in Ref. and applied to 
adiabatic quantum computing in Ref. [I3 ■ 

The A^-element search problem consists of finding a 
single marked element in a disordered A"-element list. 
The search problem is associated with an N dimensional 
Hilbert space with orthonormal basis {\k)}^^i, where the 
marked_ item corresponds to £ Following 
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where 



we consider the family of Hamiltonians 



H{s) = -il-sm{^\-s\p){pl 



N 



(1) 



(2) 



fc=i 



and s = t/T G [0,1], T being the run-time of the search. 
If the evolution is adiabatic and we start in the energy 
ground state this family of Hamiltonians takes us to 
the marked state and thus solves the search problem. 
The only relevant subspace is spanned by \tp) and We 
denote the instantaneous eigenvalues and orthonormal 
eigenvectors of H{s) restricted to the relevant subspace 
by En{s) and |£'„(s)), respectively, where n = 0, 1. We 
further define 



A(.).E.M-E,M.,/ ' + ''^-'"^-"' (3) 



N 



and 



Z{s) = \{E,{s)\EM)\ 



A useful property of Z is 



l + (iV-l)(2s-l)2 



Z{s)ds < |, 



(4) 



(5) 



for all N. 

Decoherence in the instantaneous energy eigenbasis is 
modeled by the master equation 
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-p{s) = ~tAT[H{s),p{s)] 



-BT[W{s),[W{s),p{s)]], (6) 
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where A> and B > are constants independent of N. 
Here, W{s) is assumed to be Hermitian, nondegenerate, 
and fulfill [W{s),H{s)] — 0. Furthermore, let w„(s) be 
the eigenvalues of W{s) corresponding to the eigenvectors 
\Enis)) and define r(s) = wi{s) — wo{s). 

Next, we implement the idea of local adiabatic 
search by making a monotone, sufficiently smooth 
reparametrization s € [0,1] ^ r — f{s) G [0, 1] of H{s) 
and W{s) in such a way that more time is spent near the 
minimum energy gap. In the closed case (i? = 0), it was 
shown in Ref. Q that the optimal choice 
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A2(r' 
1 



-dr\ 



'0 A2(r') 

TT N 



dr' — 



N 



arctan ( \/ N — 1) 



2 ViV" 



(7) 



yields the criterion T ^ \/N for the run-time, in 
analogy with the Grover search [T^. Applying this 
reparametrization results in the transformation s — > r 
as well as in multiplication by "^^^^ (r) of the right-hand 
side of Eq. ©. Assume that p(0) = \Eq{Q)) {Ea{Q)\ and 
let 

Y{r) = {E^{r)\p{r)\E^{r)) - {E^{r)\p{r)\E^{r)) 

= PooW - PiiW- (8) 

We now address the main objective of this paper, which 
is to determine how the probability to remain in the 
ground state depends on the run-time T and the param- 
eters A and B. The strategy is to express this probabil- 
ity, indirectly in terms of Y{r), as an integral equation. 
Thereafter, we apply appropriate estimates to obtain a 
lower bound for the probability. 

We may rewrite Eq. J^Jl as an integral equation that 
takes the form 



where 

I{r) 
u±{r') 

and 

Q{r) = 

R{t) = 
We further define 



1 - Y{r) = 4/(r), 



(9) 



'■")±^^^('-")lZ(r")r(r")dr"(10) 



^ ' dr' L 



r2(r') 



A2(r') 
1 



A(r') 



-dr'. (11) 



C = min 



rG[0,l] A(r) 



(12) 



In the case where A > 0, we wish to estimate 
|1 — F(r)|. This can be done by calculating an up- 
per bound for |/±(r)|, using Z(r)A(r) < y^{N - l)/N, 
exp[-TBQ(r)] < 1, exp{-TB[Q(r) - Q{r')]} < 1 if 
r > r' , as well as Eqs. Q and (|12|) . which result in 



|/±WI < 



ttL IN-1 

Y 



1 



Lit 
T 2 



Z(r')A(r') 



(13) 



d_ 

dr' 



B 



A{r'] 



± iA 



dr'. 



By use of Z(r)A2(r) = ^/N - l/N and Eq. lO, the in- 
tegral on the right-hand side of Eq. (|13|l can be estimated 
as 



< 







/' 
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dr' 1 




A 



Z(/)A(/) 



^2^2 + ^2 

B 

+ 52^2 +y^2 



A(r') ^ 

\{r') 



dr' 



dr' ^ ' 



dr' 



Z(r') 



,Y\r') 



dr' . (14) 



Note that we have extended the integration interval from 
[0,r] to [0,1]. Since both Z(r) and A(r) are symmetric 
around r = i, it follows that Z(r)|^A(r)| has the same 
symmetry. Moreover, A(r) is increasing on the inter- 
val [i,l]. Hence, Z(r)|^A(r)| = Z{r)lA{r) on [i,l], 
which leads to 



Z{r) 



d . , s 
-A(r) 
dr 



dr ^ 2 



< 2 



d 

/ Z(r)4-A(r)dr 
ii/2 a?' 



iV- 1 



N 



< 2. 



(15) 



To deal with the second term on the right-hand side of 
Eq. H14|l . we introduce the following condition 



Z{r) 



^T'ir) 
dr 



dr < K, 



(16) 



where i^T is a constant independent of N Since 

Br^{r) can be seen as the instantaneous strength of 

the decoherence, the condition in Eq. H16(l essentially 

states that the fluctuations in strength arc not allowed 

to grow with N. If one assumes that r(r) — 77(A(r)), 

where rj : (0, 00) (0, 00) is an increasing, sufhciently 

smooth function, it can be shown that it is sufficient that 

ri{x) < Cx'^, where C and tr > ^ are constants, to fulfill 

the condition in Eq. (|16|l . This means that the condition 

is fulfilled for the particular case where W{r) = H{r). 

By combining Eqs. Q, H13|l - and using that 

^2^2 _^ ^2 > ^2 g^jj^ 52p4(o) + ^2 > ^2^ ob^j^ijj 




(17) 
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Hence, it is a sufficient condition for local adiabaticity 
that T ^ ^/N . In conclusion, an increased degree of 
eigenbasis decoherence does not change the asymptotic 
behavior of the run-time of the adiabatic search. This 
result is independent of the explicit form of W{r) as long 
as Eq. (0 is fulfilled and [W{r),H{r)] = 0. 

In the wide-open case (A = 0), the protective effect 
of the Hamiltonian dynamics is absent and one may ex- 
pect that the asymptotic behavior depends on the explicit 
choice of W{r). To verify this point, we let W{r) be such 
that r(r) = A'^(r), o- > 1 Ha. We further put B = 1 
for convenience. Notice that the choice W{r) = H{r) 
corresponds to T{r) = A{r). We prove that in the wide- 
open case with r(r) = A'^(r), a sufficient and necessary 
condition for adiabaticity is T N"' . Note that we have 
to show that the sufficient condition is also necessary, as 
we wish to prove that the wide-open case is essentially 
different from the A case. 

To prove the sufficiency, we insert T{r) = A°'{r) into 
Eq. (HTJ, and use that A(r) > l/y/N and cr > 1, to 
obtain 

Q(r) = i A''^-\r')dr' > jj^r (18) 
Inserting Eq. ((TH|l into Eq. 1(1171) gives 

1-Y{r) < 4 f f e-TW^('-'-'-"^Z{r')Z{r")dr"dr'. 
Jo Jo 

(19) 



Finally, we use Z{r) < \/N — 1 and Eq. © to obtain 



Pooir) > 1 - 



2 T 



(20) 



Thus, a sufficient condition for adiabaticity is T ^ N"' . 

Now we show that this condition is also necessary. We 
let Yt denote the solution of Eq. ^ for a given run-time 
T. It can be proved that 0| 



Yrir) > Yo{r), 



(21) 



which means that an evolution with non-zero run-time 
remains closer to the instantaneous ground state than 
the evolution with zero run-time. Insert Eq. H21|) into 
Eq. H10() and combine with Eq. ^ to obtain 



1->t(1)>/o, 



(22) 



where 



/o = 4/ / e-^[QM-Q('^')lz(r)Z(r')ro(r')^i^'d^ 



(23) 



and 



Yo{r) = 



> - 



1 - (iV- l)(2r- 1) 
Ny/1 + {N -l){2r-iy 
{N - l){2r - 1) 



N^l + (TV - l)(2r- 1)2 



(24) 



In order for Yr(l) 1, -^o has to go to zero, since Yrir) < 
1. Hence, we have found a necessary condition for the 
system to approach adiabaticity. 

In order to express this necessary condition in terms 
of the run-time T, let us use Eq. in Eq. if^ and 
make the change of variables x = \/N — l(2r — 1) and 
y = ViV - l(2r' - 1). This yields 



h > 



N- 1 
N 



/(a,VA^), 



(25) 



where 



/(a,/3) 



a'S'iv) 



P g-a*(a:) 

,^TT^7_^, (1+2/2)3/2 



dydx (26) 



with 



$(x) 



T 



T 



2LN''-WN - 1 2N'^ arctan(VA^ - 1) 
{l + x'Y~^dx'. (27) 



Furthermore 
d 



d(3 



I{a,(3) = -2- 



i$(/3) 



P y sinh(Q;$(y)) 

—dy 



1+/32 7„ (l + y2)3/2 



+2 



/3e- 



-Q$(/3) 



(l-h/32)3/2 7, 



cosh(a<I>(a;)) 



dx 



/□ -Q*(/3) i-P -ait>(x) 
> 2-- — — / T^dx 



(l+/32)3/2 7o 1 + 

F(a,/3) > 0. 



(28) 



This expression is obtained by separating the integrals 
''I-fj" iiito + Jq" and making the change of vari- 
ables X — > —X and y ^ —y in the "/°^" integrals, as well 
as by using the inequality 

y ^ P 1 



< 



(l + j^2)(3/2) - ^/l^l + y2^ 



Vye[0,/3]. (29) 



It follows from Eq. H28I) that I{a, (3) is increasing in /3, 
which together with Eqs. (|25|l and 128() gives 



lo > -^/(a, xMT^T) > _/(a, 1) 

1 

> -^J^ F{a,(3')d(3' >0, 



(30) 



where we have assumed that N > 2. Thus, if /q ^ then 
/q F{a, P')d(3' — > necessarily. Furthermore, we have 

r-l 

F{a,(3)dl3 



d 

da ,/Q 



= -2 



(l + /32)3/2 7o l + x- 



f3e 



-Q$(/3) 



(l + /?2)3/2 7o 



$(x)e 



-dx 



-dx 



d(3 < 0. 



(31) 
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FIG. 1: Local search with W(r) = H{r) and success proba- 
bihty 0.5. The curves show logj T vs logj A*', where T is the 
run-time needed to obtain the success probabihty 0.5, and 
where A'^ is the hst length. Each curve shows the resuh for 
a given degree of decoherence u. Counted from below, the 
curves correspond to = 0, 0.1, . . . , 0.9, 0.95, 1, interpolating 
between the closed {u> — 0) and wide-open (w = 1) case. As 
seen, all curves tend to the slope i, except the uppermost 
wide-open case, which tends to the slope 1. 

It follows that /g F{a, I3)d(3 is a strictly decreasing func- 
tion in a. Hence, a necessary condition for this expression 
to go to zero is that a — > oo. For large N it follows from 
the expression for a in Eq. H27() that it is necessary for 
adiabaticity that T > A^"^. 

In Fig. ^ we supplement the above analytic results 
with numerical simulations of the dynamics of Eq. © 
with the choice W{r) — H{r) and initial condition 



Poo(O) = 1. We interpolate between the closed and 
the wide-open case, by letting A = cos(w7r/2) and 
B — sin(a)7r/2), where lo goes from to 1. Furthermore, 
we have assumed the success probability /5oo(l) = 0.5. 
These simulations confirm the predictions concerning the 
asymptotic behavior, viz., that the evolution of the local 
adiabatic quantum search stays near the instantaneous 
ground state if T ^ \/N for all cases except the wide- 
open one, where T ^ N. 

In conclusion, we have demonstrated that local adi- 
abatic search is robust to decoherence in the instanta- 
neous eigenbasis of the search Hamiltonian, as long as 
the Hamiltonian dynamics is present. Up to a condition 
on the fluctuations, this result is independent of the ex- 
plicit form of the decoherence term. This independence 
does no longer hold in absence of the Hamiltonian part, 
in which case the asymptotic behavior of the run-time 
of the local search changes. The protective effect of the 
Hamiltonian dynamics is an indication of robustness of 
quantum adiabatic search, which may be of importance 
in physical implementations of a working search scheme 
that outperforms any known classical search algorithm. 

An interesting extension would be to apply the present 
analysis to adiabatic algorithms designed to solve other 
problems, such as, e.g^ the NP-complete problems 3-SAT 
y and exact cover ]2, E| • Although analytical re- 
sults may not be achievable for these problems, numerical 
investigations could reveal whether or not the protective 
effect of the Hamiltonian dynamics is present. More- 
over, one might consider whether there occurs a tran- 
sition from the seemingly polynomial behavior found in 
y , to an exponential time-complexity, as the strength of 
decoherence increases. 

We wish to thank Patrik Thunstrom for useful com- 
ments on the manuscript. 
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